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Data from e"*'e~ annihilation into hadrons at centre-of-mass energies between 91 GeV and 209 GeV are used 
to study tiie four-jet rate as a function of the Durham algorithm's resolution parameter ycut- The four-jet rate 
is compared to next-to- leading order calculations that include the resummation of large logarithms. The strong 
coupling measured from the four-jet rate is 

a^Mz) = 0.1208 ± 0.0006(stat.) ± 0.0021(ea;p.) ± 0.0019(/wd.) ± 0.0024(scoZe) 
in agreement with the world average. 



1. Introduction 

The annihilation of an electron and a positron 
into hadrons allows a precise test of Quantum 
Chromodynamics (QCD). Multijet rates are pre- 
dicted in perturbation theory as functions of the 
jet-resolution parameter, with one free param- 
eter, the strong coupling as- Calculations be- 
yond leading order are made possible by theoret- 
ical developments achieved in the last few years. 
For multi-jet rates as well as numerous event 
shape distributions with perturbative expansions 
starting at 0{as), matched next-to-leading order 
(NLO) and next-to-leading logarithmic approxi- 
mations (NLLA) provide very precise description 
of the data over a wide range of the available kine- 
matic region and centre-of-mass energy [1-4]. 
In this analysis we use data collected by the 
OPAL Collaboration at LEP at 13 centre-of-mass 
energy points covering the full LEP energy range 
of 91-209 GeV. A more detailed description of the 
analysis can be found in [5]. 

2. Observable 

Jet algorithms are applied to cluster the large 
number of particles of an hadronic event into a 
small number of jets, reflecting the parton struc- 
ture of the event. For this analysis we use the 



Durham scheme [2]. The next-to-leading order 
calculation predicts the four-jet rate i?4, which 
is the fraction of four-jet events, as a function of 
as [6] . The fixed-order perturbative prediction is 
not reliable for small values of j/cut- An all-order 
resummation, given in Ref. [2], is possible for the 
Durham clustering algorithm. This NLLA calcu- 
lation is combined with the NLO-prediction using 
the modified "R matching" scheme described in 
Ref. [6]. 

3. Analysis Procedure 

The data used in this analysis were collected by 
OPAL [7] between 1995 and 2000 and correspond 
to 14.5 pb~^ of 91 GeV data, 10.4 pb~^ ofLEPl.5 
data, with centre-of-mass energies y/s—13Q GeV 
and 136 GeV and 706.4 pb"^ of LEP2 data with 
-^5=161 to 209 GeV. The highest energy points 
have a spread of 1-2 GeV and are grouped into 
the main energy points. The data at 91 GeV were 
taken during calibration runs at the Z*^ peak at 
the beginning of each year and during data taking 
in 1996-2000. The breakdown of the data sam- 
ples, mean ^/s and collected luminosity are given 
in Table 1. To study fewer energy points with 
improved statistical power, we combine the data 
above 91 GeV into three energy points by taking 
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luminosity times hadronic cross-section weighted 
averages of the samples. The LEP1.5 data sam- 
ples provide an energy point at 133 GeV, while 
the LEP2 samples give energy points at 178 GeV 
and 198 GeV corresponding to the energy ranges 
from 161 to 183 GeV and from 189 to 209 GeV. 
Samples of Monte Carlo simulated events (JET- 
SET 7.4[8], KK2f[9], PYTHIA 6. 125 [10], HER- 
WIG 6.2[11], KORALW 1.42[12] and GRC4f[13]) 
were used to correct the data for experimental 
resolution, acceptance and backgrounds. The 
Monte Carlo samples generated at ciacli energy 
point studied were processed through a full sim- 
ulation of the OPAL detector [14], and recon- 
structed in the same way as the data. In ad- 
dition, for comparisons with the corrected data, 
and when correcting for the effects of fragmen- 
tation, large samples of generator-level Monte- 
Carlo events were employed, using the parton 
shower models PYTHIA 6.158, HERWIG 6.2 and 
ARIADNE 4.11 [15]. 

3.1. Selection of Events 

The selection of events for this analysis con- 
sists of three main stages: the identification of 
hadronic event candidates, the removal of events 
with a large amount of initial state radiation 
(ISR), and the removal of four-fermion back- 
ground events (at energies above the W+W~ pro- 
duction threshold electroweak four-fermion pro- 
duction becomes a substantial background). 
The number of selected candidate events ob- 
tained after applying the selection cuts are given 
in Table 1. After all cuts, the acceptance for 
non-radiative signal events ranged from 88.5% 
at 91 GeV to 76.5% at 207 GeV. The resid- 
ual four-fermion background was negligible below 
161 GeV, and increased from 2.1% at 161 GeV to 
6.2% at 207 GeV. 

3.2. Corrections to the data 

The C!xpectc!<:l number of residual four-fermion 
background events, was subtracted from the num- 
ber of data events. The effects of detector reso- 
lution, acceptance and of residual ISR were then 
accounted for by a bin-by-bin correction proce- 
dure. The Monte Carlo ratio of the hadron level 
to the detector level was used as a correction fac- 



tor for the data. 
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Table 1 

The average centcr-of-mass energy and integrated 
luminosity for each data sample, together with 
the numbers of selected data and expected Monte 
Carlo events. The horizontal lines indicate the 
four energy intervals. 



4. Systematic Uncertainties 

Several sources of possible systematic uncer- 
tainties are studied. We investigate experimental 
uncertainties, effects due to hadronisation and the 
theoretical uncertainty, associated with missing 
higher order terms in the theoretical prediction. 

5. Results 

5.1. Four-Jet Rate Distributions 

The four-jet rates at y^=198 GeV after sub- 
traction of background and correction for detec- 
tor effects is shown in Figure 1. Superimposed we 
show the predictions by the PYTHIA, HERWIG 
and ARIADNE Monte Carlo models, which in all 
cases were tuned to OPAL data at 91 GeV centre- 
of-mass energy. In order to make a more sensitive 
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Figure 1. The figures show the four-jet rate dis- 
tribution at hadron-level as a function of the j/cut 
resolution parameter obtained with the Durham 
algorithm. 



retical prediction is then multiplied by the ratio 
of the hadron- and parton-level four-jet rates. 
A x^-value for each energy point is calculated. A 
single event can contribute to several t/cut-bins in 
a four-jet rate distribution and for this reason the 
bins are correlated and the complete covariance 
matrix was determined. The value is mini- 
mized with respect to as for each energy point 
separately. The scale parameter is set to unity. 
The systematic error was determined as described 
in Section 4. Each distribution was re-fitted and 
the difference to the default as value was taken 
as a symmetric systematic error. The result for 
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Table 2 

The mean value of ag for each energy interval, the 
statistical and experimental errors, and the error 
due to hadronization and scale uncertainties. 



comparison between data and models, the inserts 
in the upper right corner show the differences be- 
tween data and each model, divided by the com- 
bined statistical and experimental error in that 
bin. The three models are seen to describe the 
high energy data at ^/s= 136, 178 and 198 GeV 
well. However, some minor discrepancies are seen 
in the data taken at 91 GeV. 

5.2. Determination of ccg 

Our measurement of the strong coupling as is 
based on the fits of QCD predictions to the cor- 
rected four-jet rate distribution. The combined 
C'(as)-I-NLLA calculation was used. The theoret- 
ical predictions of the four-jet rate provide distri- 
butions at the parton level. In order to confront 
the theory with the hadron-level data, it is nec- 
essary to correct for hadronization effects. The 
four-jet rate was calculated at hadron and parton 
level using PYTHIA and, as a cross-check, with 
the HERWIG and ARIADNE models. The theo- 



Q^s(-\/s) and the systematic uncertainties within 
each energy range are summarized in Table 2. It 
is also of interest to combine the measurements of 
as from the four different centre-of-mass energy 
points in order to determine a single value. This 
problem has been subject of extensive study by 
the LEP QCD working group [16], and we adopt 
their procedure here. 

The set of as measurements to be combined 
are first evolved to a common scale, Q = M^o, 
assuming the validity of QCD. The measure- 
ments are then combined in a weighted mean, 
to minimize the between the combined val- 
ues and the measurements. The result is 
as(Mz) = 0.1208 ±0.0006(stat.)±0.0021(ea;p.)± 
0mi9{had.) ± 0.0024(scaZe), consistent with the 
world average value of 0.1172 ± 0.0020 [17]. The 
average value of as obtained and the results at 
each energy point are shown in Figure 2. The 
hadronization and the scale uncertainty decrease 
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Figure 2. The values for as in the various energy 
intervals. The errors show the statistieal (inner 
part) and the total error. The statistical error at 
91 GeV is very small and cannot be seen. 



with increasing energy. The hadronization scales 
as an inverse power of ^/s and the scale uncer- 
tainty varies with ag so that, for example, a 10% 
smaller value of as at higher energies leads to a 
27% smaller scale uncertainty. 

6. Summary 

In this note we present the preliminary mea- 
surements of the OPAL collaboration of the 
strong coupling from the four-jet rate at centre- 
of-mass energies between 91 and 209 GeV. The 
predictions of the PYTHIA, HERWIG and ARI- 
ADNE Monte Carlo models are found to be in 
general agreement with the measured distribu- 
tions. Some differences are noted in the 91 GeV 
data where the statistics are largest. 
From a fit of ©(qD+NLLA predictions to the 
four-jet rate, we have determined the strong cou- 
pling as- The value of as{Mz) is determined to 
be as(Mz) = 0.1208 ± 0.0038 (total error). 
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